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The topic of this letter is structure formation with topological defects. We 
first present a partially new, fully local and gauge invariant system of pertur- 
bation equations to treat microwave background and dark matter fluctuations 
induced by topological defects (or any other type of seeds). We show that this 
treatment is extremly well suited for linear numerical analysis of structure 
formation by applying it to the texture scenario. Our numerical results cover 
a larger dynamical range than previous investigations and are complementary 
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!_i ■ since we use substantially different methods. 
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Despite great effort and considerable progress, the problem of structure formation 
in the universe remains basically unsolved. Observations show that the fluctuation 
spectrum on the large scales observed by COBE should be not very far from scale 
invariant |IJ, |2|, |3| . This has been considered as great success for inflationary models 
which predict a scale invariant fluctuation spectrum. In this letter we consider an 
alternative class of models which also yield a scale invariant spectrum of Cosmic Mi- 
crowave Background (CMB) fluctuations: Models where perturbations are seeded 
by global topological defects which can form during symmetry breaking phase tran- 
sitions in the early universe [|J. To be specific, we consider texture, 7r 3 -defects which 
lead to event singularities in four dimensional spacetime ||. A common feature of 
global topological defects is the behaviour of the energy density in the scalar field 
which scales like px oc l/(at) 2 and thus represents always the same fraction of the 
total energy density of the universe. 

p T /p ~ 8nGr] 2 = 2e , (1) 

where r\ determines the symmetry breaking scale. The background spacetime is a 
Friedmann-Lemaitre universe with Q — 1. We choose conformal coordinates such 

that 

ds 2 = a 2 (—dt 2 + Sijdx l dx^) . 

Numerical analyses of CMB fluctuations from topological defects on large scales have 
been performed in || [7| . A spherically symmetric approximation is discussed in . 
Results for intermediate scales are presented in ||. All these investigations (except 
the quite rough spherically symmetric calculation ||) use linear cosmo logical pertur- 
bation theory in synchronous gauge and take into account only scalar perturbations. 
In this letter we derive a fully gauge invariant and local system of perturbation 
equations. The (non local) split into scalar, vector and tensor modes on the hyper- 
surfaces of constant time is not performed. We solve the equations numerically in 
a cold dark matter (CDM) universe with global texture. In this letter, we present 
the main results. In a longer paper we give detailed derivations of the equations and 



fully describe our numerical methods |T0[ . Since there are no spurious gauge modes 



in our initial conditions, their is no danger that these may grow in time and much 
of the difficulties to choose correct initial conditions (see e.g. [0]) are removed. 

We find that previously ignored vector and tensor fluctuations contribute ap- 
proximately 25 % to the total microwave background anisotropies (see Fig. 1). 

We calculate the microwave background anisotropies on angular scales which are 
larger than the angle which subtends the horizon scale at decoupling of matter and 
radiation, 6 > 6d- For Q = 1 and Zd ~ 1000 

9 d = l/ViTTT ~ 0.03 « 2° . (2) 

It is therefore sufficient to study the generation and evolution of microwave back- 
ground fluctuations after recombination. During this period, the photons are decou- 
pled from baryonic matter and free stream, influenced solely by cosmic gravitational 



redshift and by perturbations in the gravitational field (if the medium is not reion- 
ized). The photon distribution function which lives in seven dimensional relativistic 
phase space PqM. = {(x,p) G TJ\4\g(x)(p,p) = 0}, obeys Liouville's equation 

X g {f) = . (3) 

In a tetrad basis e^, X g is given by (see e.g. fill) 



X, HP% W>K|-) , (4) 

where u u are the connection 1-forms of (M,g) in the basis e M . 

The metric of a perturbed Friedmann universe with Q — 1 is given by ds 2 = 
g^dx^dx" with 

9fiu = a 2 ^^ +h la/ ) = a 2 ~g iiV , (5) 

where (r/^) = diag(— ,+,+,+) is the flat Minkowski metric and \h,j. v \ <C 1 is a 
small perturbation. We now use the fact that the motion of photons is conformally 
invariant. Taking into account the different affine parameters, (|3p is equivalent to 

{X s f){x,ap) = 0. (6) 

If e M is a tetrad in Minkowski space, e M = e M + {l/2)h v e v is a tetrad w.r.t the 
perturbed geometry g. For (x, e^) G Pq, thus, (x, e M p M ) G Po- We can therefore 
define the perturbation of the distribution function F by 

f(x, p%) = />, p%) + F(x, p%) . (7) 

Furthermore, we set p % = pY, with p 2 = X^=i(p*) 2 and v = ap. Liouville's equation 
for / then yields a perturbation equation for F. We choose the natural tetrad 
e /i = d fi + (i/2)h"pd v . Using (|,|) and (§) we obtain 

(<9 t + j%)F = -[H L + (A n +^)f + {H VJ - B hJ )YY]v^ , (8) 

where we have parametrized 

, . / 2A 25, \ 

^-^2S« 211^ + 211^ J' W 

with F* = 0. 

We now define 

m = 1/4 / Fv 3 dv . 

p r a A J 

4m is the fractional perturbation of the brightness l, 
i = a~ A / fv 3 dv . 



Setting i = i(T(j,x)), one finds that m corresponds to the fractional perturbation 
in the temperature, 

T(j,x)=f(l + m{j,x)) . (10) 

A more explicit derivation of (|T0| ) is given in Jl2|). Integrating ([8]) v 3 dv , we obtain 

d t m + Ydtm = H L + (A,i +^B i )j i + (H tj - B hj ) 7 y • (H) 

It is well known that the equation of motion for photons only couples to the Weyl 
part of the curvature (null geodesies are conformally invariant). The r.h.s. of fllTD is 
given by first derivatives of the metric only which could at most represent integrals 
of the Weyl tensor. To obtain a local, non integral equation, we thus rewrite ( |Ti~D in 
terms of Am. In fact, defining 

X = Am - (AH L - itf,g) - l -AB a l , 
( |Tl~j) yields for \ the equation of motion 

(d t + 7 ^)x = -Zrf&Eij - 'f-feuABy , (12) 

where €ku is the totally antisymmetric tensor in three dimensions, £7y and Bij are the 
electric and magnetic part of the Weyl tensor. The spatial indices in this equation 
are raises and lowered with Sij and therefore no care is taken in the index position. 
Double indices are summed over irrespective of their position. 

In eqn. ([E^) the contribution from the electric part of the Weyl tensor does not 
contain tensor perturbations. On the other hand, scalar perturbations do not induce 
a magnetic gravitational field. The second contribution to the source term in fll2|) 
represents a combination of vector and tensor perturbations. If vector perturbations 
are negligible, the two terms on the r.h.s of (|12"D therefore represent a split into scalar 
and tensor perturbations which is local. 

In terms of metric perturbations, the electric and magnetic part of the Weyl 
tensor are given by (see, e.g. [|Kj) 

Eij = \[A^A - H L ) - b l3 - (Afly + \iitki) ~ Hi, -Hf lH } (13) 

■&ij " 7)1 ^ilm'-'jmil itjlm "imil ) ? V / 
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with <Jij = ~(Bi,j +Bj,i ) - -5 i:j Bi - Hij and A {j = didj - (1/3)^- A . 
Since the Weyl tensor of Friedmann Lemaitre universes vanishes, the rhs of (|T2| ) 



is manifestly gauge invariant. (This is the so called Stewart lemma fl6].) Therefore 
also the variable \ is gauge invariant. 



The general solution to ([12|) is given by 



X (t, x, 7) = / S T {t', x + (t' - t) 7 , j)dt' + X (ti, x + (ti - t) 7 , 7) , (15) 



where St is the source term given on the rhs of dl2|) . 

The electric and magnetic part of the Weyl tensor are determined by the per- 
turbations in the energy momentum tensor via Einstein's equations. We assume 
that the source for the geometric perturbations is given by the scalar field and dark 
matter. The contributions from radiation may be neglected. Furthermore, vector 
perturbations of the dark matter (which decay quickly) are neglected. The diver- 
gence of Eij is then determined by (see, |L3| or JTJJ]) 

&Eij = -MGpDMl'Di - BnGidiSToo + 3(-)ST oi ) + 12ttG3%- , (16) 

CI 

where 

T .. = T .. _ (a 2 /3)^-T/ = r% exture) = 0, 4 0„ -(1/3)<%(V0) 2 , 

dT 0j = dT 0j = 4>(f> d , 

6T 00 = 5Tl ) t exture) = 1 -m 2 + (V<f ) ) 2 ), 

and Dj is a gauge invariant perturbation variable for the density gradient (see [|T^, 
14, l0|1). For scalar perturbations Dj = djD. The evolution equation for the dark 



matter density perturbation is given by 

D + (-)D - AnGa 2 p DM D = 8ttG0 2 . (17) 



The equation for B^ is more involved. A somewhat cumbersome derivation [ |10| 
yields 



B tJ + 3(^)B y - AB l3 = SnGS^ 



:(B) 



with S {j = ei m (i5T i,j) m —(a/a)ei m (iTj)i, m . 

Here (i...j) denotes symmetrization in the indices i and j. 

To these equations we have to add the evolution equation of the scalar field, 

<j> + 2(d/o)0 - A0 = a 2 A0(0 2 - rj 2 ) . (19) 

We have solved the closed hyperbolic system flT2"| , [IB], [Tj], [18], and |TUp numerically 
on a 192 3 grid for different initial conditions on a NEC-SX3 computer at the Centro 
Svizzero di Calcolo Scientifico (CSCS). The numerical methods employed and the 



different tests of our programs are described in |pl| . Here we just want to present 
the main results. 

Since on subhorizon scales gauge dependent and gauge invariant variables do not 
differ substantially we can interpret the variables D and \ by 

D = (Sp/p) = 5 and X = A(<JT/f ) . 

Using fast Fourier transforms we calculate the spectrum P(k) = \S(k)\ 2 and 5T/T, 
which we then expand in spherical harmonics 

(5T/f)(t ,x,~f) =J2 a Ux)Y lm (7) ■ (20) 

I 111. 

As usual we assume that the average over different observer positions coincides with 
the ensemble average and determine 

c - = wrwIt laUx)l2 - (21) 

We have performed 10 simulations on a 192 3 grid with about 100 different observer 
positions for each simulation. The average harmonic amplitudes with ler variance 
are shown in Fig. 2. The low order multipoles depend strongly on the random 
initial conditions (cosmic variance), like in the spherically symmetric simulation ||. 
From Fig. 2 it is clear that the texture scenario is compatible with a scale invariant 
spectrum. The main difference of the currently favored scenarios with inflation 
induced perturbations lies in the distribution of fluctuations which is non Gaussian 
in models with topological defects. The pixel distribution for AT/T in the sky is 
negatively skewed. Our quadrupole apmlitude is given by 

Q = (0.53±0.16)e 

To reproduce the COBE amplitude Qcobe = (0.6 ± 0.1)10~ 5 [|3|, we have to nor- 
malize our spectrum by choosing the phase transition scale r\ 

e = 47rGr] 2 = (1.1 ± 0.5)10" 5 . (22) 

This value is comparable with the value of e obtained in || [7]]. 

The power spectrum of dark matter density fluctuations is shown in Fig. 3. To 
be compatible with observations (a 8 « 1), we have to introduce a somewhat high 
bias factor of 

6^4±2 . 

(The bias factor takes into account that the observed clustering of light does not 
necessarily coincide with the clustering of the underlining dark matter distribution.) 
Observations and simulations of nonlinear clustering of dark matter and baryons 
|T7| hint that a bias factor b=2 - 2.5 might be reasonable. 

In this letter we have presented a closed system of cosmological perturbation 
equations which are not plagued by gauge modes and which is well suited for nu- 
merical analysis. Our numerical results are consistent with previous investigations 



H |7|, |] , indicating that the texture scenario of structure formation yields somewhat 
too much power on small scales. 
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Figure Captions 



Fig. 1 

The amplitude of the electric and magnetic source terms to the photon equation if 
motion are shown as a function of wavenumber k in arbitrary scale. On very large 
scales the magnetic part contributes about 1/4 decaying to roughly 1/10 on small 



scales (E = ± £»(^-) 2 , B = \ E^jik&B 
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Fig. 2 

The harmonic amplitudes 1(1 + 1)q are shown with la error bars. The slight rise at 
large / is due to the finite size of the texture core. No smoothing is applied. 



Fig. 3 

The power spectrum of the CDM fluctuations induced by texture. (The vertical 
scale is arbitrary.) 
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